Valence bond solids for SU(n) spin chains: 
exact models, spinon confinement, and the Haldane gap 
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To begin witii, we introduce several exact models for SU(3) spin chains: (1) a translationally 
invariant parent Hamiltonian involving four-site interactions for the trimer chain, with a three-fold 
degenerate ground state. We provide numerical evidence that the elementary excitations of this 
model transform under representation 3 of SU(3) if the original spins of the model transform under 
rep. 3. (2) a family of parent Hamiltonians for valence bond solids of SU(3) chains with spin reps. 
6, 10, and 8 on each lattice site. We argue that of these three models, only the latter two exhibit 
spinon confinement and hence a Haldane gap in the excitation spectrum. We generalize some of our 
models to SU(n). Finally, we use the emerging rules for the construction of VBS states to argue 
that models of antiferromagnetic chains of SU(n) spins in general possess a Haldane gap if the spins 
transform under a representation corresponding to a Young tableau consisting of a number of boxes 
A which is divisible by n. If A and n have no common divisor, the spin chain will support deconfined 
spinous and not exhibit a Haldane gap. If A and n have a common divisor different from n, it will 
depend on the specifics of the model including the range of the interaction. 

PACS numbers: 75.10.Jm, 75.10.Pq, 75.10.Dg, 32.80.Pj 



I. INTRODUCTION 

Quantum spin chains have been a most rewarding sub- 
ject of study almost since the early days of quantum me- 
chanics, beginning with the invention of the Bethe ansatz 
in 1931 [1] as a method to solve the 5 = ^ Heisenberg 
chain with nearest-neighbor interactions. The method 
led to the discovery of the Yang-Baxter equation in 
1967 [2, 3], and provides the foundation for the field of 
integrable models [4]. Faddeev and Takhtajan [5] dis- 
covered in 1981 that the elementary excitations (now 
called spinous) of the spin- 1/2 Heisenberg chain carry 
spin 1/2 while the Hilbert space is spanned by spin flips, 
which carry spin 1. The fractional quantization of spin 
in spin chains is conceptually similar to the fractional 
quantization of charge in quantized Hall liquids [6, 7]. In 
1982, Haldane [8, 9] identified the 0(3) nonlinear sigma 
model as the effective low-energy field theory of SU(2) 
spin chains, and argued that chains with integer spin pos- 
sess a gap in the excitation spectrum, while a topological 
term renders half-integer spin chains gapless [10, 11]. 

The general methods — the Bethe ansatz method and 
the use of effective field theories including bosoniza- 
tion [12, 13] — are complemented by a number of ex- 
actly solvable models, most prominently among them 
the Majumdar-Ghosh (MG) Hamiltonian for the S = \ 
dimer chain [14], the AKLT model as a paradigm of the 
gapped S" = 1 chain [15, 16], and the Haldane-Shastry 
model (HSM) [17-20]. The HSM is by far the most so- 
phisticated among these three, as it is not only solvable 
for the ground state, but fully integrable due to its Yan- 
gian symmetry [20]. The wave functions for the ground 
state and single-spinon excitations are of a simple Jas- 
trow form, elevating the conceptual similarity to quan- 
tized Hall states to a formal equivalence. Another unique 



feature of the HSM is that the spinous are free in the 
sense that they only interact through their half-Fermi 
statistics [21-25], which renders the model an ideal start- 
ing point for any perturbative description of spin systems 
in terms of interacting spinous [26]. The HSM has been 
generahzed from SU(2) to SU(n) [27-32]. 

For the MG and the AKLT model, only the ground 
states are known exactly. Nonetheless, these models have 
amply contributed to our understanding of many aspects 
of spin chains, each of them through the specific concepts 
captured in its ground state [33-51]. The models are 
specific to SU(2) spin chains. We will review both models 
below. 

In the past, the motivation to study SU(n) spin sys- 
tems with n > 2 has been mainly formal. The Bethe 
ansatz method has been generalized to multiple compo- 
nent systems by Sutherland [52], yielding the so-called 
nested Bethe ansatz. In particular, this has led to a 
deeper understanding of quantum integrability and the 
applicability of the Bethe ansatz [53]. Furthermore, the 
nested Bethe ansatz was used to study the spectrum of 
the SU(n) HSM [22, 27]. It has also been applied to 
SU(2) spin chains with orbital degeneracy at the SU(4) 
symmetric point [54, 55]. Most recently, Damerau and 
Kliimper obtained highly accurate numerical results for 
the thermodynamic properties of the SU(4) spin-orbital 
model [56]. SU(n) Heisenberg models have been studied 
recently by Kawashima and Tanabe [57] with quantum 
Monte Carlo, and by Paramckanti and Marston [58] using 
variational wave functions. 

The effective field theory description of SU(2) spin 
chains by Haldane yielding the distinction between gap- 
less half-integer spin chains with deconfined spinous and 
gapped integer spin chains with confined spinous cannot 
be directly generalized to SU(n), as there is no direct 



equivalent of the CP^ representation used in Haldane's 
analysis. The critical behavior of SU(n) spin chains, how- 
ever, has been analyzed by Affleck in the framework of 
effective field theories [59, 60]. 

An experimental realization of an SU(3) spin system, 
and in particular an antiferromagnetic SU(3) spin chain, 
however, might be possible in an optical lattice of ul- 
tracold atoms in the not-too-distant future. The "spin" 
in these systems would of course not relate to the ro- 
tation group of our physical space, but rather relate to 
SU(3) rotations in an internal space spanned by three de- 
generate "colors" the atom may assume, subject to the 
requirement that the number of atoms of each color is 
conserved. A possible way to realize such a system ex- 
perimentally is described in Appendix A. Moreover, it 
has been suggested recently that an SU(3) trimer state 
might be realized approximately in a spin tetrahedron 
chain [61, 62]. 

Motivated by both this prospect as well as the math- 
ematical challenges inherent to the problem, we propose 
several exact models for SU(3) spin chains in this article. 
The models are similar in spirit to the MG or the AKLT 
model for SU(2), and consist of parent Hamiltonians and 
their exact ground states. There is no reason to expect 
any of these models to be integrable, and none of the 
excited states are known exactly. We generalize several 
of our models to SU(7i), and use the emerging rules to 
investigate and motivate which SU(n) spin chains exhibit 
spinon confinement and a Haldanc gap. 

The article is organized as follows. Following a brief re- 
view of the MG model in Sec. II, we introduce the trimer 
model for SU(3) spin chains in Sec. III. This model con- 
sists of a translationally invariant Hamiltonian involv- 
ing four-site interactions, with a three-fold degenerate 
ground state, in which triples of neighboring sites form 
SU(3) singlets (or trimers). In Sec. IV, we review the rep- 
resentations of SU(3), which we use to verify the trimer 
model in Sec. V. In this section, we further provide nu- 
merical evidence that the elementary excitations of this 
model transform under representation 3 of SU(3) if the 
original spins of the model transform under rep. 3. We 
proceed by introducing Schwingcr bosons in Sec. VI and 
a review of the AKLT model in Sec. VII. In Sec. VIII, 
we formulate a family of parent Hamiltonians for valence 
bond sohds of SU(3) chains with spin reps. 6, 10, and 8 
on each lattice site, and proof their validity. We argue 
that only the rep. 10 and the rep. 8 model, which are 
in a wider sense generalizations of the AKLT model to 
SU(3), exhibit spinon confinement and hence a Haldane- 
type gap in the excitation spectrum. In Sec. IX, wc gen- 
eralize three of our models from SU(3) to SU(n). In 
Sec. X, we use the rules emerging from the numerous 
VBS models we studied to investigate which models of 
SU(n) spin chains in general exhibit spinon confinement 
and a Haldane gap. In this context, we first review a rig- 
orous theorem due to Affieck and Lieb [63] in Sec. X A. In 
Sec. XB, we argue that the spinous in SU(?7) spin chains 
with spins transforming under reps, with Young tableaux 



consisting of a number of boxes A which is divisible by 
n are always confined. In Sec. XC, we construct several 
specific examples to argue that if A and n have a common 
divisor different from n, the model will be confining only 
if the interactions are sufficiently long ranged. Specifi- 
cally, the models we study suggest that if q is the largest 
common divisor of A and n, the model will exhibit spinon 
confinement only if the interactions extends at least to 
the ^-th neighbor on the chain. If A and n have no com- 
mon divisor, the spinous will be free and chain will not 
exhibit a Haldanc gap. We briefly summarize the dif- 
ferent categories of models in Sec. XD, and present a 
counter-example to the general rules in Sec. XE. We 
conclude with a brief summary of the results obtained in 
this article in Sec. XL 

A brief and concise account of the SU(3) VBS models 
we elaborate here has been given previously [64] . 

II. THE MAJUMDAR-GHOSH MODEL 

Majumdar and Ghosh [14] noticed in 1969 that on a 
linear spin S = ^ chain with an even number of sites, the 
two valence bond solid or dimer states 

ki?) = n (4t4u-cI,4it)io> = 

(i odd) 

{[ o — o o — o o — o o ) "even" 
(1) 
I o o — o o — o o — o ) "odd" 

where the product runs over all even sites i for one state 
and over all odd sites for the other, are exact zero-energy 
ground states [65] of the parent Hamiltonian 

Hmg = + ls,S,+2 + , (2) 

where 

^^ = 1 ^tf^TT'C^r'^ (3) 

r,T'=Ta 

and cr = {ax, a-y,az) is the vector consisting of the three 
Pauli matrices. 

The proof is exceedingly simple. We rewrite 

HmG^\Y^^' H^^{S., + S,+ , + S^+2Y (4) 
i 

Clearly, any state in which the total spin of three neigh- 
boring spins is ^ will be annihilated by Hi. (The total 
spin can only be ^ or |, as i (g) | (g) | = | ffi i ® |.) In 
the dimer states above, this is always the case as two of 
the three neighboring spins are in a singlet configuration, 
and (g) ^ = i . Graphically, we may express this as 

H,\c^ o)=H,\o c^)^0. (5) 
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As Hi is positive definite, the two zero-energy eigenstatcs 
of Hmg are also ground states. 

Is the Majumdar-Ghosh or dimer state in the univer- 
sality class generic to one-dimensional spin-i liquids, and 
hence a useful paradigm to understand, say, the nearest- 
neighbor Hcisenbcrg chain? The answer is clearly no, as 
the dimer states (1) violate translational symmetry mod- 
ulo translations by two lattice spacings, while the generic 
liquid is invariant. 

Nonetheless, the dimer chain shares some important 
properties of this generic liquid. First, the spinon 
excitations — here domain walls between "even" and 
"odd" ground states — are deconfined. (To construct ap- 
proximate eigenstates of Hmgi we take momentum su- 
perpositions of localized domain walls.) Second, there 
are (modulo the overall two-fold degeneracy) only M + 1 
orbitals available for an individual spinon if 2M spins are 
condensed into dimers or valence bond singlets. This is 
to say, if there are only a few spinous in a long chain, the 
number of orbitals available to them is roughly half the 
number of sites. This can easily be seen graphically: 



even 



odd 



If we start with an even ground state on the left, the 
spinon to its right must occupy an even lattice site and 
vice versa. The resulting state counting is precisely 
what one finds in the Haldanc-Shastry model, where 
it is directly linked to the half-Fermi statistics of the 
spinous [21]. 

The dimer chain is further meaningful as a piece of a 
general paradigm. The two degenerate dimer states (1) 
can be combined into an S* = 1 chain, the AKLT chain, 
which serves as a generic paradigm for S = 1 chains 
which exhibit the Haldane gap [8, 9, 34], and provides the 
intellectual background for several of the exact models we 
introduce further below. Before doing so, however, we 
will now introduce the trimer model, which constitutes 
an SU(3) analog of the MG model. 



states on the chain are given by 



(m) \ _ 



trimer / 



= < 



) = 



^ trimcr / ' 



i'l'L.) , (6) 



T^trimcr 



Introducing operators c,]g. which create a fermion of color 
(T (cr = b, r, g) at lattice site i, the trimer states can be 
written as 



T-^t rimer 



n ( sign(^)4^4^, 



(-3^ integer j 7r(b,r,g) 



(7) 

where /i ~ 1,2,3 labels the three degenerate ground 
states, and i runs over the lattice sites subject to the 
constraint that is integer. The sum extends over all 
six permutations tt of the three colors b, r, and g, i.e., 



sign(7r)4„c|_^i^4_^2, 

(Q,;3,7)=ir(b,r,g) 



"ibS+lr'-i+2g SrS+lgW+2 



^1 



!b + 4g4+lb'^', 



.-clJ 



ib'"i+lgS+2r '--igS-|-lrS-|-2b 



igS+lbS+2r 
4r'^i-|-lb''l+2g- 



(8) 



The SU(3) generators at each lattice site i are in anal- 
ogy to (3) defined as 



J, 



a= 1 



, . . . , 



(9) 



CT.cr'— b,r.g 



where the A"^ arc the GcU-Mann matrices (see App. B) 
The operators (9) satisfy the commutation relations 



[Jf , J}] = <5,, r^Vf , a,6,c=l. 



(10) 



(we use the Einstein summation convention) with /''^'^ 
the structure constants of SU(3) (see App. B). Wc fur- 
ther introduce the total SU(3) spin of v neighboring sites 
i,. . . A + V - I, 



III. THE TRIMER MODEL 



A. The Hamiltonian and its ground states 



r(^) 



(11) 



Consider a chain with N lattice sites, where has to 
be divisible by three, and periodic boundary conditions 
(PBCs). On each lattice site we place an SU(3) spin 
which transforms under the fundamental representation 
3, i.e., the spin can take the values (or colors) blue (b), 
red (r), or green (g). The trimer states are obtained 
by requiring the spins on each three neighboring sites to 
form an SU(3) singlet, which we call a trimer and sketch 
it by o — o — o . The three linearly independent trimer 



where Ji is the eight-dimensional vector formed by its 
components (9). The parent Hamiltonian for the trimer 
states (7) is given by 



40 
~9 



(12) 



The JiJj terms appear complicated in terms of Gell- 
Mann matrices, but are rather simply when written out 
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FIG. 1: Tensor product of two S = \ spins with Young 
tableaux and weight diagrams of the occurring SU(2) rep- 
resentations. S~ is the diagonal generator. 
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FIG. 2; Tensor product of three S — ^ spins with Young 
tableaux. For SU(n) with n > 2, the tableau with three 
boxes on top of each other exists as well. 

using the operator P^ , which permutes the SU(n) spins 
(here n = 3) on sites i and j: 



J iJ 1 



1 



Pi 



(13) 



To verify the trimer Hamiltonian (12), as well as for 
the valence bond solid (VBS) models we propose below, 
we will need a few higher-dimensional representations of 
SU(3). We review these in the following section. 



IV. REPRESENTATIONS OF SU(3) 
A. Young tableaux and representations of SU(2) 

Let us begin with a review of Young tableaux and the 
representations of SU(2). The group SU(2) has three 
generators 5", a — 1,2,3, which obey the algebra 



(14) 



where repeated indices are summed over and e°'"^ is the 
totally antisymmetric tensor. The representations of 
SU(2) are classified by the spin S', which takes integer 
or half-integer values. The fundamental representation 
of SU(2) has spin S = ^, it contains the two states |t) 
and II). Higher-dimensional representations can be con- 
structed as tensor products of fundamental representa- 
tions, which is conveniently accomplished using Young 
tableaux (see e.g. [66]). These tableaux are constructed 
as follows (see Figs. 1 and 2 for examples). For each of 
the N spins, draw a box numbered consecutively from 
left to right. The representations of SU(2) are obtained 
by putting the boxes together such that the numbers as- 
signed to them increase in each row from left to right and 





FIG. 3: (Color online) a) Weight diagram of the fundamen- 
tal SU(3) representation 3 = (1,0). b) Weight diagram of 
the complex conjugate representation 3 = (0, 1). J'* and J* 
denote the diagonal generators, f/^, and the raising 
operators. 



in each column from top to bottom. Each tableau indi- 
cates symmetrization over all boxes in the same row, and 
antisymmetrization over all boxes in the same column. 
This implies that we cannot have more than two boxes 
on top of each other. If nt denotes the number of boxes 
in the ith row, the spin is given by 5 = — K2). 

To be more explicit, let us consider the tensor product 
^ <8) ^ eg" I depicted in Fig. 2 in detail. We start with the 



state If, I) = ITTT), and hence find 



> 



ITiT)+IJTT)). (15) 



The two states with S = S^- 



must be orthogonal to 



(15). A convenient choice of basis is 



1 2 ' 2 ' 



1 

71 



(16) 



where uj = exp(i^) . The tableaux tell us primarily that 
two such basis states exist, not what a convenient choice 
of orthonormal basis states may be. 

The irreducible representations of SU(2) can be clas- 
sified through the eigenvalues of the Casimir operator 
given by the square of the total spin S'^. The special fea- 
ture of 5^^ is that it commutes with all generators S°' and 
is hence by Schur's lemma [67] proportional to the iden- 
tity for any finite-dimensional irreducible representation. 
The eigenvalues are given by 

= '^SU(2) = + 



B. Representation theory of SU(3) 

The group SU(3) has eight generators J°, a — 1, 
which obey the algebra 



(17) 
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[I] ® [2] ® [3] 



mm 
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[ij ffitmi 1 8 8 

3 6 

FIG. 4; Tensor product 3 (g) 3 (g) 3 with Young tableaux. 



/i2 columns boxes 

FIG. 5: Dynkin coordinates {fii,fi2) for a given Young 
tableau. The columns containing three boxes represent ad- 
ditional SU(3) singlet factors, which yield equivalent repre- 
sentations and hence leave the Dynkin coordinates (/^i,/i2) 
unchanged. 



where the structure constants /"^"^ are given in App. B. 
For SU(3) we have two diagonal generators, usually cho- 
sen to be and J^, and six generators which define 
the ladder operators /± = ± iJ^, = ± ij'^, 
and = J"* ± iJ^, respectively. An explicit realization 
of (17) is, for example, given by the J°'s as expressed 
in terms of Gell-Mann matrices in (9). This realization 
defines the fundamental representation 3 of SU(3) illus- 
trated in Fig. 3a. It is three-dimensional, and we have 
chosen to label the basis states by the colors blue (b), 
red (r), and green (g). The weight diagram depicted in 
Fig. 3a instructs us about the eigenvalues of the diagonal 
generators as well as the actions of the ladder operators 
on the basis states. 

All other representations of SU(3) can be constructed 
by taking tensor products of reps. 3, which is again 
most conveniently accomplished using Young tableaux 
(see Fig. 4 for an example). The antisymmctrization 
over all boxes in the same column implies that we can- 
not have more than three boxes on top of each other. 
Each tableaux stands for an irreducible representation 
of SU(3), which can be uniquely labeled by their high- 
est weight or Dynkin coordinates (/Xi,/i2) [67, 68] (see 
Fig. 5). For example, the fundamental representation 
3 has Dynkin coordinates (1,0). Note that all columns 
containing three boxes are superfluous, as the antisym- 
mctrization of three colors yields only one state. In par- 
ticular, the SU(3) singlet has Dynkin coordinates (0,0). 
In general, the dimension of a representation (/ii,^) is 
given by ^(/ii -I- l){n2 + 1)(mi + M2 + 2). The labeling 
using bold numbers refers to the dimensions of the repre- 
sentations alone. Although this labeling is not unique, it 
will mostly be sufficient for our purposes. A representa- 
tion m and its conjugated counterpart m arc related to 
each other by interchange of their Dynkin coordinates. 
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FIG. 6: Weight diagram of the representation 6 = (2, 0). The 
weight diagram of the conjugate representation 6 = (0, 2) is 
obtained by reflection at the origin [67]. 



C. Examples of representations of SU(3) 

We now consider some specific representations of SU(3) 
in detail. As starting point we use the fundamental rep- 
resentation 3 spanned by the states |b), |r), and |g). The 
second three-dimensional representation 3 is obtained by 
antisymmetrically coupling two reps. 3. The Dynkin co- 
ordinates of the rep. 3 are (0,1), i.e., the reps. 3 and 3 
arc complex conjugate of each other. An explicit basis of 
the rep. 3 is given by the colors yellow (y), cyan (c), and 
magenta (m), 

1 



|y) = -^(|rg)-|gr)). 



1 



|c) = -=(|gb)-|bg)). 



^/2 
1 



(18) 



(|br)-|rb)). 



The weight diagram is shown in Fig. 3.b. The generators 
arc given by (9) with A° replaced by — (A°)*, where * 
denotes complex conjugation of the matrix elements [68]. 
In particular, we find |y) = — |c), [/+ |c) = — |m), and 

V+\y)=-H- . 

The six-dimensional representation 6 has Dynkin coor- 
dinates (2,0), and can hence be constructed by symmet- 
rically couphng two reps. 3. The basis states of the rep. 
6 are shown in Fig. 6. The conjugate representation 6 
can be constructed by symmetrically coupling two reps. 
3. 

Let us now consider the tensor product 3 (X) 3 = 1 © 8. 
The weight diagram of the so-called adjoint representa- 
tion 8 = (1,1) is shown in Fig. 7. The states can be 
constructed starting from the highest weight state |bni), 
yielding /" jbm) = |rm), U~ |bm) = — jbc), V" |bni) = 
jgm) — |by), and so on. This procedure yields two linearly 
independent states with = = 0. The representa- 
tion 8 can also be obtained by coupling of the reps. 6 
and 3, as can be seen from the Young tableaux in Fig. 4. 
On a more abstract level, the adjoint representation is 
the representation we obtain if we consider the genera- 
tors J° themselves basis vectors. In the weight diagram 
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75(|by)"kc» 

gc> ^(|by) + |rc>-2|gm)) 



FIG. 7: Weight diagram of the adjoint representation 8 = 
(1, 1). The state with = J** = is doubly degenerate [67]. 
Note that two reps. 8 can be constructed by combining three 
fundamental reps. 3 (colors), just as two reps. | can be con- 
structed by combining three SU(2) spins (cf. (16)). The states 
in the diagram span a basis for one of these representations. 
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*7s(|bgg) + |gbg) + |ggb)) 



FIG. 8: Weight diagram of the representation 10 — (3,0). 
The weight diagram of the conjugate representation 10 — 
(0, 3) is obtained by reflection at the origin [67] . 



shown in Fig. 7, the generators and correspond to 
the two states at the origin, whereas the ladder opera- 
tors , , and correspond to the states at the six 
surrounding points. In the notation of Fig. 7, the singlet 
orthogonal to 8 is given by --i=(|by) + jrc) + |gm)). 

The weight diagrams of four other representations rel- 
evant to our purposes below are shown in Figs. 8 to 10. 

It is known that the physical properties of SU(2) spin 
chains crucially depend on whether on the lattice sites 
are integer or half-integer spins. A similar distinction 
can be made for SU(3) chains, as elaborated in Sec. X. 
The distinction integer or half-integer spin for SU(2) is 
replaced by a distinction between three families of irre- 
ducible representations of SU(3): cither the number of 
boxes in the Young tableau is divisible by three without 
remainder {e.g. 1, 8, 10, 27), with remainder one {e.g. 
3, 6, 15, 15'), or with remainder two {e.g. 3, 6, 15, 15'). 

While SU(2) has only one Casimir operator, SU(3) has 
two. The quadratic Casimir operator is defined as 



J2 = ^ JV", 



(19) 
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FIG. 9: Weight diagram of the representations 15 = (2, 1) 
and 15' = (4,0). 

27 



FIG. 10: Weight diagram of the self-conjugate representation 
27 = 

erate [67]. 



(2,2). The state with = J* = is three-fold degen- 



As commutes with all generators J" it is proportional 
to the identity for any finite-dimensional irreducible rep- 
resentation. The eigenvalue in a representation with 
Dynkin coordinates (^i,/i2) is [67] 



M2 



3^1 



(20) 

We have chosen the normalization in (19) according to 
the convention 



(adjoint representation) 



where the J°'s are the generators of the representation. 



n2 

which yields C|jj^3^(l,l) = 3 for the representation 8. 
Note that the quadratic Casimir operator cannot be used 
to distinguish between a representation and its conjugate. 
This distinction would require the cubic Casimir opera- 
tor [67] , which we will not need for any of the models we 
propose below. 



V. THE TRIMER MODEL (CONTINUED) 

B. Verification of the model 

We will now proceed with the verification of the trimer 
Hamiltonian (12). Since the spins on the individual sites 
transform under the fundamental representation 3, the 
SU(3) content of four sites is 

3® 3(g) 3® 3 = 3-3 e 2-6 e 3-15 e 15', (21) 

i.e., we obtain representations 3, 6, and two non- 
equivalent 15-dimensional representations with Dynkin 
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FIG. 11: Couplings used in the numerical studies to create 
(a) the locahzed rep. 3 trial state and (b) the localized rep. 3 
trial state. 



coordinates (2, 1) and (4, 0), respectively. All these repre- 
sentations can be distinguished by their eigenvalues of the 

quadratic Casimir operator, which is given by (■^i''^) if 
the four spins reside on the four neighboring lattice sites 
I, . . . , i + 3. 

For the trimer states (6), the situation simplifies as we 
only have the two possibilities 



o — o o — o 



= 1 (g) 3 = 3, 
= 3 (g) 3 = 3 



6, 



which implies that the total SU(3) spin on four neigh- 
boring sites can only transform under representations 3 
or 6. The eigenvalues of the quadratic Casimir operator 
for these representations arc 4/3 and 10/3, respectively. 
The auxiliary operators 



H,, = 



r(4) 



hence annihilate the trimer states for all values of i, while 
they yield positive eigenvalues for 15 or 15', i.e., all other 
states. Summing Hi over all lattice sites i yields (12). 
We have numerically confirmed by exact diagonahzation 
of (12) for chains with = 9 and 12 lattice sites that 
the three states (7) are the only ground states. 

Note that the representation content of five neighbor- 
ing sites in the trimer chains is just the conjugate of the 
above, as 

c^-o — o = 1 (g) 3 = 3, 

o o — o — o o =3(g)l(g)3 = 3e6. 

Since the quadratic Casimirs of conjugate represen- 
tations have identical eigenvalues, C|uj3j(^i,/i2) = 
^su(3) (^^2, Ml); can construct another parent Hamil- 
tonian for the trimer states (7) by simply replacing jf''^ 

with J^P in (12). This Hamiltonian will have a differ- 
ent spectrum. In comparison to the four-site interaction 
Hamiltonian (12), however, it is more complicated while 
bearing no advantages. We will not consider it further. 



C. Elementary excitations 

Let us now turn to the low-lying excitations of (12). 
In analogy with the MG model, it is evident that the 



mom 


-Btot 




% 


over- 


[2tt/N] 


exact 


trial 


off 


lap 





2.9735 


4.5860 


54.2 


0.9221 


1, 12 


6.0345 


10.2804 


70.4 


0.5845 


2, 11 


9.0164 


17.2991 


91.9 


0.0 


3, 10 


6.6863 


13.1536 


96.7 


0.0 


4, 9 


3.0896 


5.0529 


63.5 


0.8864 


5, 8 


4.8744 


7.5033 


53.9 


0.8625 


6, 7 


8.5618 


16.6841 


94.9 


0.1095 



TABLE I: Energies of the rep. 3 trial states (23) in comparison 
to the exact excitation energies of the trimer model (12) and 
their overlaps for an SU(3) spin chain with A'^ = 13 sites. 



SU(3) spinon or "coloron" excitations correspond to do- 
main walls between the degenerate ground states. For 
the trimer model, however, there arc two different kinds 
of domain walls, as illustrated by: 



3 

) o ( 

3 

o — o 



(23) 
(24) 



The first domain wall (23) connects ground state /i to 
the left to ground state /i -I- 1 to the right, where ji is 
defined modulo 3 (see (7)), and consists of an individ- 
ual SU(3) spin, which transforms under representation 
3. The second domain wall (24) connects ground state 
/i with ground state /i 4- 2. It consists of two antisym- 
metrically coupled spins on two neighboring sites, and 
hence transforms under representation 3. As we take 
momentum superpositions of the localized domain walls 
illustrated above, we expect one of them, but not both, to 
constitute an approximate eigenstate of the trimer model. 
The reason we do not expect both of them to yield a 
valid excitation is that they can decay into each other, 
i.e., if the rep. 3 excitation is valid the rep. 3 domain 
wall would decay into two rep. 3 excitations, and vice 
versa. The question which of the two excitations is the 
valid one, i.e., whether the elementary excitations trans- 
form under 3 or 3 under SU(3) rotations, can be resolved 
through numerical studies. We will discuss the results of 
these studies now. 

The rep. 3 and the rep. 3 trial states require chains 
with A = 3 • integer + 1 and A = 3 • integer + 2 sites, re- 
spectively; we chose A^ = 13 and A^ = 14 for our numeri- 
cal studies. To create the localized domain walls (23) and 
(24), we numerically diagonalizcd auxiliary Hamiltoni- 
ans with appropriate couplings, as illustrated in Fig. 11. 
From these localized excitations, we constructed momen- 
tum eigenstates by superposition, and compared them to 
the exact eigenstates of our model Hamiltonian (12) for 
chains with the same number of sites. The results are 
shown in Tab. I and Fig. 12 for the rep. 3 trial state, and 
in Tab. II and Fig. 13 for the rep. 3 trial state. 



7 




exact energy state for N=1 4 ■ 
trial state • 






momenta [2 pi/N] 



momenta [2 pi/N] 



FIG. 12: Dispersion of the rep. 3 trial states (23) in compar- 
ison to the exact excitation energies of (12) for a chain with 
A'^ = 13. The lines are a guide to the eye. 



FIG. 13: Dispersion of the rep. 3 trial states (24) in compar- 
ison to the exact excitation energies of (12) for a chain with 
N = 14. The lines are a guide to the eye. 



mom 


-Etot 




% 


over- 


[2tt/N] 


exact 


trial 


off 


lap 





2.1013 


2.3077 


9.8 


0.9953 


1, 13 


4.3677 


4.8683 


11.5 


0.9864 


2, 12 


7.7322 


8.7072 


12.6 


0.9716 


3, 11 


6.8964 


7.7858 


12.9 


0.9696 


4, 10 


3.2244 


3.5415 


9.8 


0.9934 


5, 9 


2.2494 


2.4690 


9.7 


0.9950 


6, 8 


5.4903 


6.1016 


11.1 


0.9827 


7 


7.4965 


8.5714 


14.3 


0.9562 



TABLE II: Energies of the rep. 3 trial states (24) in compar- 
ison to the exact excitation energies of the trimer model (12) 
and their overlaps for an SU(3) spin chain with A'' = 14 sites. 



reason is simply that domain walls connect one ground 
state with another, without introducing costly correla- 
tions in the region between the domain walls. In the case 
of the MG model and the trimer model introduced here, 
however, there is still an energy gap associated with the 
creation of each coloron, which is simply the energy cost 
associated with the domain wall. 

In most of the remainder of this article, we will intro- 
duce a family of exactly soluble valence bond models for 
SU(3) chains of various spin representations of the SU(3) 
spins at each lattice site. To formulate these models, we 
will first review Schwinger bosons for both SU(2) and 
SU(3) and the AKLT model. 

VI. SCHWINGER BOSONS 



The numerical results clearly indicate that the rep. 3 
trial states (24) are valid approximations to the elemen- 
tary excitations of the trimer chain, while the rep. 3 trial 
states (23) are not. We deduce that the elementary exci- 
tations of the trimer chain (12) transform under 3, that 
is, under the representation conjugated to the original 
SU(3) spins localized at the sites of the chain. Using the 
language of colors, one may say that if a basis for the 
original spins is spanned by blue, red, and green, a basis 
for the excitations is spanned by the complementary col- 
ors yellow, cyan, and magenta. This result appears to be 
a general feature of SU(3) spin chains, as it was recently 
shown explicitly to hold for the Haldane-Shastry model 
as well [30-32]. 

Note that the elementary excitations of the trimer 
chain are deconfined, meaning that the energy of two 
localized representation 3 domain walls or colorons (24) 
does not depend on the distance between them. The 



Schwinger bosons [69, 70] constitute a way to formulate 
spin-S" representations of an SU(2) algebra. The spin 
operators 

S'^ + iSy = S+ = a)b, 
S^-iSy = S- ^ b'fa, (25) 
5^ = i(ata-5t6), 

are given in terms of boson creation and annihilation op- 
erators which obey the usual commutation relations 

[a,at] = [6,5t] =1, 
^ ^ ^ ^ (26) 
[a, 6] = [a,6t] = [a\b] = [a^b^ =0. 

It is readily verified with (26) that S^, Sy , and satisfy 
(14). The spin quantum number S is given by half the 
number of bosons, 

2S^a^a + b^b, (27) 
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and the usual spin states (simultaneous eigenstates of S'^ 
and S^') are given by 

\S,m)- )J „ j J O). (28) 

In particular, the spin-i states are given by 

|T)=4|0)=at|0), U)=4|0)=&t|0), (29) 

i.e., and act just like the fermion creation opera- 
tors cj and cj in this case. The difference shows up only 
when two (or more) creation operators act on the same 
site or orbital. The fermion operators create an antisym- 
metric or singlet configuration (in accordance with the 
Pauli principle), 



|0,0)=4cj|0). 



(30) 



while the Schwinger bosons create a totally symmetric or 
triplet (or higher spin if we create more than two bosons) 
configuration. 



|f,f) = ^(at)2|o; 

|f,0) = at6t|0), 
|1,-1) = 7f(^^)'|0) 



(31) 



The generalization to SU(n) proceeds without incident. 
We content ourselves here by writing the formalism out 
explicitly for SU(3). In analogy to (25), we write the 
SU(3) spin operators (9) 







= /+ 


= b'fr, 


Jl 




= I- 










^ l{b'fb-rfr), 


J4 


+ iJ^ 


- v+ 




J4 




= V- 


= g^b, 


J6 




= u+ 


= r^g, 


J6 




= u- 


= g^r, 






J8 


= 5^(&^& + rtr-2.gt.g), 



(32) 



in terms of the boson annihilation and creation operators 
b,b^ (blue), r,r^ (red), a.nd g,g^ (green) satisfying 



[Ml 



r, r' 



(33) 



while all other commutators vanish. Again, it is readily 
verified with (33) that the operators J° satisfy (17). The 
basis states spanning the fundamental representation 3 
may in analogy to (29) be written using either fermion 
or boson creation operators: 



(34) 



|b) 




0) 


= 6t 





|r) 




0) 


= rt 





|g) 


= 4 


0) 


-g' 


|0 



We write this abbreviated 

3 = (1,0) = □ = ctjO) = at|0; 



(35) 



The fermion operators can be used to combine spins 
transforming under the fundamental representation 3 an- 
tisymmetrically, and hence to construct the representa- 
tions 



3 
1 



(0,1) 
(0,0) 



B 



44 1 o; 



(36) 



— c^clc^g I ) 



The Schwinger bosons, by contrast, combine fundamental 
representations 3 symmetrically, and hence yield repre- 
sentations labeled by Young tableaux in which the boxes 
are arranged in a horizontal row, like 



6 = (2,0) = m = «^/3^|0), 
10 = (3,0) = nJJ = at/3V|0) 



(37) 



15' = (4,0) = 



/3t7t<5t|o), 



where a, /?, 7, . . . G {&, r, g}. Unfortunately, it is not pos- 
sible to construct representations like 
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(1,1) 



by simply taking products of anti-commuting or commut- 
ing creation or annihilation operators. 



VII. THE AKLT MODEL 

Using the SU(2) Schwinger bosons introduced in the 
previous section, we may rewrite the Majumdar-Ghosh 
states (1) as 



(38) 



This formulation was used by Affleck, Kennedy, Lieb, 
and Tasaki [15, 16] to propose a family of states for higher 
spin representations of SU(2). In particular, they showed 
that the valance bond solid (VBS) state 



IV'aklt) = n - |0) = 

i 

= vI/^gK'&I • *mgK'&^] |0) 



o — o 
o o — o 



o — o o — o o 
o)-^ o — o o — o 

T 

projection onto spin 5=1 



(39) 
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is the exact zero-energy ground state of the spin-1 ex- 
tended Heisenberg Hamiltonian 

HaKLT = (s^S^+l + \{S,S,+ iY + ^^j (40) 



with periodic boundary conditions. Each term in the 
sum (40) projects onto the subspace in which the total 
spin of a pair of neighboring sites is 5 = 2. The Hamil- 
tonian (40) thereby lifts all states except (39) to positive 
energies. The VBS state (39) is a generic paradigm as it 
shares all the symmetries, but in particular the Haldane 
spin gap [8, 9, 34], of the spin-1 Heisenberg chain. It even 
offers a particularly simple understanding of this gap, or 
of the linear confinement potential between spinous re- 
sponsible for it, as illustrated by the cartoon: 



VIII. SU(3) VALENCE BOND SOLIDS 

To begin with, we use SU(3) Schwinger bosons intro- 
duced in Sec. VI to rewrite the trimer states (7) as 

4±.) = n ( E «ig^('^)«I/3l+i7l+2)|0) 

N (".'3,7)= 
(-TT- integer j TT{b,r,g) 

= *^[fet,rt,5t] |o), (43) 

where, as in (7), = 1,2,3 labels the three degenerate 
ground states, i runs over the lattice sites subject to the 
constraint that is integer, and the sum extends over 
all six permutations tt of the three colors b, r, and g. 
This formulation can be used directly to construct VBSs 
for SU(3) spin chains with spins transforming under rep- 
resentations 6 and 10 on each site. 



o — o o — o (p O O O O (|) o — o o — o 

o — o o — o o — o o — o o — o o — o o — o 



energy cost oc distance 



Our understanding [71, 72] of the connection between 
the confinement force and the Haldane gap is that the 
confinement effectively imposes an oscillator potential for 
the relative motion of the spinous. We then interpret the 
zero-point energy of this oscillator as the Haldane gap in 
the excitation spectrum. 

The AKLT state can also be written as a matrix prod- 
uct [42, 44, 45]. We first rewrite the valence bonds 



and then use the outer product to combine the two vec- 
tors at each site into a matrix 




10) 



|1,0), a/2|1,-1) 
-V2|l,l), -|1,0), 



(41) 



Assuming PBCs, (39) may then be written as the trace 
of the matrix product 



A. The representation 6 VBS 

We obtain a representations 6 VBS from two trimer 
states by projecting the tensor product of two funda- 
mental representations 3 onto the symmetric subspace, 
i.e., onto the 6 in the decomposition 3 (g) 3 = 3 © 6. 
Graphically, this is illustrated as follows: 

one site 



T 

projection onto rep. 6 = (2, 0) 



(44) 



This construction yields three linearly independent 6 
VBS states, as there are three ways to choose two dif- 
ferent trimer states out of a total of three. These three 
VBS states are readily written out using (43), 



6 VBS / 



lo; 



(45) 



for /i = 1, 2, or 3. If we pick four neighboring sites on 
a chain with any of these states, the total SU(3) spin of 
those may contain the representations 



= 3 (g) 3 = 3 



6 



or the representations 



o — o o — o 



= 3 



3 = 2 • 3 ® 6 ® 15, 



i.e., the total spin transforms under 3, 6, or 15 = (1,2), 
all of which are contained in the product 



I^aklt) =tr(]^A/, 



(42) 



In the following section, we will propose several exact 
models of VBSs for SU(3). 



6(8)6«)6®6 = 

3-3©6-6©7-T5©3-T5'©3-21 

© 8 • 24 © 6 • 42 © 45 © 6 • 60 © 3 • 63 (46) 

and hence possible for a representation 6 spin chain 
in general. The corresponding Casimirs are given by 
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n2 

'-'SU(3) 



(0,1) 



3' ^SU(3)(^'0) — 3°' ^'^'^ '^SU(3)(-'^' ~ 3^- 

This leads us to propose the parent Hamiltonian 



6VBS 



TV 

E 



with 



H,, = 



(47) 



16 



(48) 

Note that the operators Jf , a = 1, . . . , 8, are now given 
by 6 X 6 matrices, as the Gell-Mann matrices only provide 
the generators (9) of the fundamental representation 3. 
Since the representations 3, 6, and 15 possess the small- 
est Casimirs in the expansion (46), Hi and hence also 
^6VBS E^re positive semi-definite (i.e., have only non- 
negative eigenvalues). The three linearly independent 
states (45) are zero-energy eigenstates of (47). 

To verify that these are the only ground states, we have 
numerically diagonalized (47) for = 6 and iV = 9 sites. 
For A'^ = 9, we find zero-energy ground states at momenta 
k = 0,3, and 6 (in units of ^ with the lattice constant 
set to unity). Since the dimension of the Hilbert space 
required the use of a LANCZOS algorithm, we cannot 
be certain that there are no further ground states. We 
therefore diagonalized (47) for = 6 as well, where we 
were able to obtain the full spectrum. We obtained five 
zero energy ground states, two at momentum k = and 
one each at fc = 2,3,4. One of the ground states at 
k = and the fc = 2, 4 ground states constitute the space 
of momentum eigenstates obtained by Fourier transform 
of the space spanned by the three 6 VBS states (45). 
The remaining two states at fc = 0, 3 are the momentum 
eigenstates formed by superposition of the state 




i = 6 1 2 3 4 5 6 1 ^^^^ 
and the same translated by one lattice spacing. It is 
readily seen that these two states are likewise zero energy 
eigenstates of (47) for = 6 sites. The crucial difference, 
however, is that the 6 VBS states (45) remain zero-energy 
eigenstates of (47) for all A^'s divisible by three, while 
the equivalent of (49) for larger A^ do not. We hence 
attribute these two additional ground states for A^ = 6 
to the finite size, and conclude that the three states (45) 
are the only zero-energy ground states of (47) for general 
A^'s divisible by three. 

Excitations of the 6 VBS model are given by domain 
walls between two of the ground states (45). As in the 
trimer model, two distinct types of domain walls exist, 
which transform according to representations 3 and 3: 




vl>^ .vl> 



It is not clear which excitation has the lower energy, and 
it appears likely that both of them are stable against 
decay. Let us first look at the rep. 3 excitation. The 
four-site Hamiltonian (48) annihilates the state for all i's 
except the four sites in the dashed box in (50), which 
contains the representations 

3(8)3(g)3®3®3 

= 6- 3©5-6©6-T5©T5'©2-24e42 

i.e., the representations is' = (0,4), 24 = (3,1) twice, 
and 42 = (2, 3) with Casimirs ^ and respectively, 
in addition to representations annihilated by Hi . For the 
rep. 3 excitation sketched on the right in (50), there are 
two sets of four neighboring sites not annihilated by Hi 
as indicated by the dashed and the dotted box. Each set 
contains the representations 



3 = 3-3©3-6©215 



24 



i.e., only the rep. 24 in addition to representations anni- 
hilated by Hi. For our parent Hamiltonian (47), it hence 
may well be that the rep. 3 anti-coloron has the lower en- 
ergy, but it is all but clear that the rep. 3 has sufficiently 
higher energy to decay. For general representation 6 spin 
chains, it may depend on the specifics of the model which 
excitation is lower in energy and whether the conjugate 
excitation decays or not. 

Since the excitations of the rep. 6 VBS chain are 
merely domain walls between different ground states, 
there is no confinement between them. We expect the 
generic antiferromagnetic rep. 6 chain to be gapless, even 
though the model we proposed here has a gap associated 
with the energy cost of creating a domain wall. 



B. The representation 10 VBS 

Let us now turn to the 10 VBS chain, which is a 
direct generalization of the AKLT chain to SU(3). By 
combining the three different trimer states (43) for /i = 
1,2, and 3 symmetrically, 

1^10 VBs) = [b\ r\ .gt] .^-2 [^t^ ^t] .^3 [^f ,4 ^t] | q ) 
= n( E sign(7r)at4^^^t A|o), (51) 



(50) 



■I (q,/3,7) = 

we automatically project out the rep. 10 in the decom- 
position 3(g)3(g)3 = 1©2-8©10 generated on each 
lattice site by the three trimer chains. This construction 
yields a unique state, as illustrated: 

one site 



(52) 




projection onto 10 = (3, 0) 
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In order to construct a parent Hamiltonian, note first 
that the total spin on two (neighboring) sites of a rep. 
10 chain is given by 



10 (g) 10 = 10 © 27 ® 28 e 35. 



(53) 



On the other hand, the total spin of two neighboring sites 
for the 10 VBS state can contain only the representations 

3®3®3®3 = 2-le4-8e 10 0100 27, (54) 

as can be seen easily from the dashed box in the cartoon 
above. (Note that this result is independent of how many 
sites we include in the dashed box.) After the projection 
onto rep. 10 on each lattice site, we find that only reps. 
To = (0, 3) and 27 =■ (2, 2) occur for the total spin of 
two neighboring sites for the 10 VBS state. With the 
Casimirs C|u(3)(0,3) = 6 and C|u(3)(2,2) = 8 we obtain 
the parent Hamiltonian 



N 

10 VBS = ^ 

4=1 



(55) 



the operators Jf, a ~ 1,...,8, are now 10 x 10 ma- 
trices, and we have used = 6. i?iovBS is positive 
semi-definite and annihilates the 10 VBS state (51). We 
assume that (51) is the only ground state of (55). 

The Hamiltonian (55) provides the equivalent of the 
AKLT model [15, 16], whose unique ground state is con- 
structed from dimer states by projection onto spin 1, for 
SU(3) spin chains. Note that as in the case of SU(2), it is 
sufficient to consider linear and quadratic powers of the 
total spin of only two neighboring sites. This is a general 
feature of the corresponding SU(n) models, as we will 
elaborate in the following section. 

Since the 10 VBS state (51) is unique, we cannot have 
domain walls connecting different ground states. We 
hence expect the coloron and anti-coloron excitations to 
be confined in pairs, as illustrated below. The state be- 
tween the excitations is no longer annihilated by (55), 
as there are pairs of neighboring sites containing higher- 
dimensional representations, as indicated by the dotted 
box below. As the number of such pairs increases lin- 
early with the distance between the excitation, the con- 
finement potential depends linearly on this distance. 



coloron 
3 



anti-coloron 
3 




energy cost oc distance 



(56) 



In principle, it would also be possible to create three col- 
orons (or three anti-colorons) rather than a coloron-anti- 
coloron pair, but as all three excitations would feel strong 
confinement forces, we expect the coloron-anti-coloron 
pair to constitute the dominant low energy excitation. 
The confinement force between the pair induces a linear 



oscillator potential for the relative motion of the con- 
stituents. The zero-point energy of this oscillator gives 
rise to a Haldane-type energy gap (see [71, 72] for a sim- 
ilar discussion in the two- leg Heisenberg ladder), which 
is independent of the model specifics. We expect this 
gap to be a generic feature of rep. 10 spin chains with 
short-range antiferromagnetic interactions. 



C. The representation 8 VBS 

To construct a representation 8 VBS state, consider 
first a chain with alternating representations 3 and 3 on 
neighboring sites, which we combine into singlets. This 
can be done in two ways, yielding the two states 



c^O c^O c^—O 
3 3 



and 



O — o O — o O 



We then combine a 3-3 state with an identical one shifted 
by one lattice spacing. This yields representations 3 
3 = 1 8 at each site. The 8 VBS state is obtained by 
projecting onto representation 8. Corresponding to the 
two 3-3 states illustrated above, we obtain two linearly 
independent 8 VBS states, and 'f^, which may be 
visualized as 



one site 



o — O c^O 



and 



o — o lo — o ; o — o 



projection onto 8 = (1, 1) 



(57) 



These states transform into each other under space re- 
flection or color conjugation (interchange of 3 and 3). 

It is convenient to formulate the corresponding state 
vectors as a matrix product. Taking (b,r,g) and (y,c,m) 
as bases for the reps. 3 and 3, respectively, the singlet 
bonds in above can be written 



(|b).Jy).+i 



= (|b).,|r)^.,|g). 



/ |y>.+i 

VHm, 



We are hence led to consider matrices composed of the 
outer product of these vectors on each lattice site. 




|b).,|r),,|g), 



In the case of the AKLT model reviewed above, the 
Schwinger bosons take care of the projection automat- 
ically, and we can simply assemble these matrices into 
a product state. For the 8 VBS, however, we need to 
enforce the projection explicitly. This is most elegantly 
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accomplished using the Gell-Mann matrices, yielding the 
projected matrix 

= -^Anr(A'^Mi®^). (58) 

a=l 

Here we have simply used the fact that the eight Gell- 



Mann matrices, supplemented by the unitary matrix, 
constitute a complete basis for the space of all complex 
3x3 matrices. By omitting the unit matrix in the ex- 
pansion (58), we effectively project out the singlet state. 
Written out explicitly, we obtain 



/l|by>. 



re) 



|bc), 
|bm). 



4lby). 



ky>. 



|rm). 



gm; 



ilby). 



|gy>i 
|gc). 



(59) 



|gm)j / 



Assuming PBCs, the 8 VBS state '^^ (illustrated in on 
the left (57)) is hence given by the trace of the matrix 
product 



IV'svBs) =tr n^^' 



(60) 



To obtain the state (illustrated on the right in (57)) 
we simply have to transpose the matrices in the product. 



I V's VBS 



(61) 



Let us now formulate a parent Hamiltonian for these 
states. If we consider two lattice sites on an SU(3) chain 
with a representation 8 on each lattice site in general, we 
find the full SU(3) content 



8(8)8 = 102-80 10 e 10 ©27 



(62) 



with 10 (3,0), 10 = (0,3), and 27 = (2,2). On the 
other hand, for the 8 VBS states only the representations 
3 (8) 3 = 1 ® 8 can occur for the total spin of two neigh- 
boring sites, as the two sites always contain one singlet 
(see dashed box in (57) on the right above). With the 
Casimirs C|u(3^(0,0) = and C|u(3^(l, 1) = 3 for repre- 
sentations 1 and 8, respectively, we construct the parent 
Hamiltonian 



N 



8 VBS 



(63) 



where the operators Jf , a = 1, . . . , 8, are now 8x8 ma- 
trices, and we have used the Casimir J ^ = 3 on each 
site. -ffsvBS is positive semi-definite, and annihilates 
the states and . We have numerically verified for 
chains with = 3, 4, 5, and 6 lattice sites that ^E''^ and 
are the only ground states of (63). 
Naively, one might assume the 8 VBS model to support 
deconfined spinous or colorons, which correspond to do- 
main walls between the two ground states and 4*^. A 



closer look at the domain walls, however, shows that this 
is highly unlikely, as each domain wall is a bound state 
of either two anti-colorons or two colorons, as illustrated 
below. 

anti-colorons colorons 
3 3 3 3 

o — O o — O • O-o O — o O-o • o — O 
o — O • O — o O^ O — o • cy-O o-O 



(64) 



There is no reason to assume that the domain wall de- 
picted above as two anti-colorons in fact corresponds to a 
single coloron, as it appears to be the case for the trimer 
chain. There we created a domain wall corresponding 
to a single coloron by removing one of the rep. 3 spins 
from a trimer, leaving the remaining rep. 3 spins coupled 
antisymmetrically as in the ground state. If we were to 
combine the two reps. 3 into a rep. 3 in (64), we would 
not reproduce a correlation present in the ground state, 
but enforce a new correlation. The correct interpretation 
of the domain wall between and 4"^ is hence that of a 
bound state between two linearly confined anti-colorons. 
The origin of the confining potential is illustrated below. 



anti-coloron 
3 



anti-coloron 
3 



^>—0 o — O c^O o — O • O — o 
o — O • O — o O — o O — o O^ 



energy cost oc distance 



(65) 



As in the 10 VBS, the confinement induces a linear 
oscillator potential for the relative motion of the anti- 
colorons. The zero-point energy of this oscillator cor- 
responds to a Haldane-type gap in the spectrum. The 
ground state wave function of the oscillator is symmet- 
ric, and hence corresponds to a symmetric combination 
of 3 (8 3, i.e., rep. 6. The antisymmetric combination 
3 corresponds to the first excited state of the oscillator. 
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rep. 10 □ 
rep. 27 o 
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k [2tt/N] 

FIG. 14: Spectrum of the 8 VBS Hamiltonian (63) for Af = 8 
sites obtained by exact diagonalization. (The hnes are merely 
guides to the eye.) The "magnon" excitation transforming un- 
der rep. 8 of SU(3) has the lowest energy, followed by a singlet 
excitation, as expected from the discussion in the text. The 
well defined modes at low energies provide strong evidence of 
coloron-anti-coloron bound states as compared to deconfined 
domain walls, and hence support our conclusion that the 8 
VBS exhibits a Haldane gap due to spinon confinement. 

which we expect to cost more than twice the energy of 
the symmetric state [72]. This statement holds for the 
pair of colorons in (64) as well. 

The domain walls, however, are not the only low en- 
ergy excitations. In either of the ground states, we can 
create coloron-anti-coloron bound states, which make no 
reference to the other ground state, as illustrated below. 

anti-coloron coloron 
3 3 

o — O o — O o — O o — O o — O o — O 

c^O • O-o • c^KD 

energy cost oc distance (^QQ) 

The oscillator model tells us again that the "symmetric" 
combination of 3 (X) 3, i.e., rep. 8, has the lowest energy, 
which we expect to be comparable, if not identical, to 
the energy required to create each of the domain walls 
above. The singlet 1 will have an energy comparable to 
that of a domain wall transforming under cither 3 or 3. 
In any event, we expect the 8 VBS model to display a 
Haldane gap due to coloron confinement. 

The excitation spectrum of (63) for a chain with = 8 
sites and PBCs is shown in Fig. 14. The spectrum shows 
that the lowest excitation transforms under rep. 8, as 
expected from (66), with a singlet and then another 
rep 8 following at slightly higher energies. It is tempt- 
ing to interpret those three levels as the lowest levels of 



the coloron-anti-coloron oscillator (66), but then there 
should be another singlet at a comparable spacing above. 
The fact that the spacings between these excitations are 
significantly smaller than the energy of the first exited 
state, however, would be consistent with such an inter- 
pretation, as the spinous in VBS models always have a 
local energy cost associated with their creation, which is 
specific to these models and not related to the universal 
Haldane gap stemming from confinement forces. 

Most importantly, the spectrum provides strong evi- 
dence in favor of our assumption that the domain walls 
are not elementary excitations, but bound states of ei- 
ther two colorons or two anti-colorons, and hence that 
the lowest energy excitations of finite chains are coloron- 
anti-coloron bound states as illustrated in (66). The as- 
sumption is crucial for our conclusion that the model 
exhibits a Haldane gap. If the low energy sector of the 
model was determined by two deconfined domain walls, 
we would see a continuum of states in the spectrum, sim- 
ilar to the spectrum seen in spin S = ^ chains of SU(2). 
The well defined low-energy modes in Fig. 14, however, 
look much more like the spinon-spinon bound state exci- 
tations seen in 5 = 1 chains or two-leg S = Heisenberg 
ladders. In particular, if we assume that the individual 
domain walls transform under reps. 6 and 6, we expect 
excitations transforming under the representations con- 
tained in 6 (g) 6 = 1 8 © 27 to be approximately degen- 
erate. Fig. 14 shows clearly that such a multiplet is not 
present a the lowest energies. 

IX. SU(n) MODELS 

In this section, we generalize three of the models pro- 
posed for SU(3) spin chains, the trimer model, the sym- 
metric representation 10 VBS, and the matrix product 
state 8 VBS to the case of SU(7i) spin chains. 

A. The n-mer model 

Consider an SU(n) spin chain with N sites, where 
is a multiple of n, with a spin transforming according 
to the fundamental representation n of SU(7i) at each 
lattice site, 

n = (1,0,..., 0) = □ = 4|0), (67) 

where tj denotes a "flavor", ct G {fi, . . . , f„}, and cj. cre- 
ates a fcrmion of flavor a. 

The SU(n) generators at site i are in analogy to (3) 
and (9) defined as 

Jt^l E 4.V:a'C,^', a = 1, . . . , n^-l, (68) 

where the V"' denote the n? — 1 SU(n) Gell-Mann ma- 
trices [73]. The generators are normalized through the 
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eigenvalue the quadratic Casimir operator takes in the 
adjoint representation, = C|u(„)(li 0, . . . , 0, 1) = n. 

To determine the eigenvalues of the quadratic Casimir 
for general representations of SU(n), a significant amount 
of representation theory is required [74] . We content our- 
selves here by providing the formulas up to n = 6 in 
App. C. 

In analogy to the trimcr states (7), we construct the 
n-mer states of an SU(n) spin chain by combining sets of 
n neighboring spins into a singlet, 

n 

V'iier) = n ( E «ign(^) n 1 0), 

i {a\,...,(7n)— — 1 

integer) 7r{fi , . . .,f„) 

(69) 

where /i = 1, . . . , n labels the n degenerate ground states, 
and i runs over the lattice sites subject to the constraint 
that is integer. The sum extends over all n\ permu- 
tations TT of the n flavors fi, . . . , f„. 

In order to identify a parent Hamiltonian, consider the 
total SU(n) spin on n -I- 1 neighboring sites for the n-mer 
states. Following the rules of combining representations 
labeled by Young tableaux (see e.g. [67, 68]), it is not 
difficult to see that the total spin will only contain rep- 
resentations given by tableaux with n -I- 1 boxes and two 
columns, i.e., tableaux of the form 




V rows 



with \ < V < -^^Y^- The eigenvalues of the quadratic 
Casimir operator for these representations are 

fn{i^) = ^ (^1^2^ - 1) - 2n{u - If - l) . (70) 

An educated guess for a parent Hamiltonian for the ri- 
mer chain hence appears to be 

N 

Htrial = ^^Hi (71) 

with 




where [ J denotes the floor function, i.e., \x\ is the largest 
integer I < x, and we use the notation introduced in (11). 

This construction yields the MG model [14] for SU(2), 
the trimer model (12) for SU(3), and a valid parent 
Hamiltonian for the four degenerate 4-mer states for 
SU(4). For ri > 5, however, the decomposition of 
the tensor product ri'*'"+^) contains irreducible repre- 
sentations corresponding to Young tableaux with more 



than two columns, whose Casimirs are equal or smaller 
than a number of Casimirs included in the list /n(i'), 
1/ = 1,2, ... , [^^Y^J. If the Casimir of such an "unde- 
sired" representation not included in the list is smaller 
than an odd number of Casimirs included in the list, we 
obtain negative eigenvalues for Hi, and it is not a pri- 
ori clear any more that the Hamiltonian (71) is positive 
semi-definite. An obvious cure to this problem is to write 

N 

H„-mci = ^ Hf, (73) 

i=l 

with Hi as in (72). This does, however, not cure po- 
tential problems arising from undesired representations 
which share the eigenvalues of the Casimir with one of 
the representations from the list, as it happens to be the 
case for n = 5. The Hamiltonian (73) likewise annihilates 
these representations, giving rise to a remote possibility 
that the n-mer states (69) are not the only ground states 
of (73). The potential relevance of these problems has to 
be investigated for each n separately. 

B. The representation (n, 0, . . . , 0) VBS 

As a generalization of the AKLT model for SU(2) and 
the 10 VBS model for SU(3) discussed above, we now 
consider a VBS for an SU(n) chain of spins transforming 
under the symmetric representation 

(n,0,...,0) = II II = blbl...bljO), 
n boxes 

where each a S {fi, . . . , f„}, is an SU(n) Schwinger 
boson. The VBS state is obtained by combining n n-mer 
states (69), one for each = 1, . . . ,7i, in that the total 
spin on each lattice site is projected onto the symmetric 
representation (n, 0, . . . , 0). This yields 

n 

|^(»,o,...,o)VBs) - n( E sign(7r) ]^6j_i_^^^^J|0). 

i ((Ti cr,i)= K=l 

7r(fi,...,f„) 

(74) 

Let us now construct a parent Hamiltonian for the 
symmetric VBS (74). The total SU(n) spin of two neigh- 
boring sites of a representation (n, 0, . . . , 0) spin chain in 
general contains all the representations corresponding to 
Young tableaux with 2n boxes and at most two rows, i.e., 
all tableaux of the form 




n — v columns 2u boxes 

The eigenvalues of the quadratic Casimir operator for 
these representations are given by 

g„(i/) = C2u(„)(2i/, n-j^, 0, . . . , 0) = 2n2 - An + v{v + 1). 

(75) 
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On the other hand, the total SU(n) spin of two neigh- 
boring sites of the representation (n, 0, . . . , 0) VBS (74) 
has to be contained in the product 



B 



□ ® □ 



(76) 



As we project the spin on each lattice site onto the rep- 
resentation (n, 0, . . . , 0), only two these representations 
remain: 



and 



re- 



The eigenvalues of the quadratic Casimir operator are 
given by 5„(0) = 2n{n - 2) and 5„(1) = 
spectively. Hence, using J1 = n{n — 1) 
parent Hamiltonian 



we obtain the 



(n,0,...,0) VBS 

N 



^((J,J,+i) + {2n-l)J,J,+i+n{n-\) 



(77) 



Since .gri(O) > for n > 2 and gn{i^) is a strictly in- 
creasing function of v, the Hamiltonian (77) is positive 
semi-definite. For n = 2, we recover the AKLT model 
(40); for n = 3, we recover the 10 VBS model (55). 



C. An example of a matrix product state 

In principle, a matrix product VBS can be formu- 
lated on all SU(n) chains with spins transforming under 
the symmetric combination of any representation m and 
its conjugate representation m. Unless the rep. m is 
self-conjugate, we obtain two inequivalent states, which 
transform into each other under space reflection. The 
construction of these is analogous to the 8 VBS intro- 
duced above, and likewise best illustrated as 

one site 



m — m 
-m m- 



m- 
-rfi 



m — m 
-rn m- 



and 



projection onto the symmetric 
combination in m ® m 



(78) 



The thick lines indicate that wc combine pairs of neigh- 
boring representations m and rn into singlets. On 
each lattice site, we project onto the symmetric com- 
bination of m and m, as indicated. By "symmetric 
combination" we mean that if representations m and 
m of SU(n) have Dynkin coordinates (/^i,/i2, . ■ . ,/in-i) 
and (^„_i, //„_2, ■ . ■ , ^i), respectively, we combine them 
into the representation with Dynkin coordinates {^i + 
^J'n-ll ^J'2^^J'7l-2, ■ ■ ■ iMn-i+A*!)- In other words, we align 
the columns of both tableaux horizontally, and hence ob- 
tain a tableau with twice the width, without ever adding 



a single box vertically to a column of the tableaux we 
started with. The states (78) we obtain are translation- 
ally invariant and we expect the parent Hamiltonians to 
involve nearest-neighbor interactions only. 

In this subsection, we will formulate the simplest 
SU(n) model of this general family. We take m to be 
the representation formed by antisymmetrically combin- 
ing K < ^ fundamental representations. 



m = [k] = (0, 



. ,0,1,0,. ..,0) = 

I 

K-th entry 



K boxes 



which implies that we consider a model with the represen- 
tation corresponding to a Young tableaux with a column 
with n — K boxes to the left of a column with k boxes at 
each lattice site: 



The construction of the parent Hamiltonian is simi- 
lar to the n-mer model above. The total spin on two 
neighboring lattice sites can only assume representations 
contained in m ®rn, i.e., representations corresponding 
to tableaux of the form 



with Q < u < K. The eigenvalues of the quadratic Casimir 
operator for these representations are 

K{v) = u{n-u + 1). (79) 

The obvious proposal for a parent Hamiltonian is hence 



N 



i/=0 



H = Y,m, i^. = n -K{i^)), (80) 



where \x\ denotes again the floor function. This Hamil- 
tonian singles out the matrix product state (78) as unique 
ground states for n < 5, but suffers from the same short- 
comings as (71) with (72) for n > 6. 



X. SPINON CONFINEMENT AND THE 
HALDANE GAP 

A. The Affleck— Lieb theorem 

For the generic SU(2) spin chain, Haldane [8, 9] identi- 
fied the 0(3) nonlinear sigma model as the effective low 
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energy field theory of SU(2) spin chains, and argued that 
chains with integer spin possess a gap in the excitation 
spectrum, while a topological term renders half-integer 
spin chains gapless [10, 11]. The exact models for SU(2) 
spin chains we reviewed above, the MG and the AKLT 
chain, serve a paradigms to illustrate the general princi- 
ple. Unfortunately, the effective field theory description 
of Haldanc yielding the distinction between gapless half 
integer spin chains with deconfined spinous and gapped 
integer spin chains with confined spinous cannot be di- 
rectly generalized to SU(7i) chains, as there is no direct 
equivalent of the CP^ representation used in Haldane's 
analysis. 

Nonetheless, there is a rigorous and rather general re- 
sult for antiferromagnetic chains of SU(7i) spins trans- 
forming under a representation corresponding to a Young 
tableau with a number of boxes not divisible by n: Af- 
fleck and Licb [63] showed that if the ground state is 
non-dcgencratc, and the Hamiltonian consists nearest- 
neighbor interactions only, than the gap in the excita- 
tion spectrum vanishes as 1/7V (where N is the number 
of sites) in the thermodynamic limit. This result is fully 
consistent with the picture suggested by the models in- 
troduced above. Like the MG model, the trimer model 
and the representation 6 VBS have degenerate ground 
states and interactions which extend beyond the nearest 
neighbor, which implies that the theorem is not directly 
applicable. 

On physical grounds, however, the statement that a 
given model is gapless {i.e., the excitation gap vanishes 
in the thermodynamic limit) implies that the spinous are 
deconfined. The reason is simply that if there was a con- 
finement force between them, the zero-point energy asso- 
ciated with the quantum mechanical oscillator of the rel- 
ative motion between the spinous would inevitably yield 
an energy gap. The MG, the trimer and the 6 VBS model 
constitute pedagogically valuable paradigms of decon- 
fined spinons. Since the excitations in these models are 
literally domain walls between different ground states, no 
long-range forces can exist between them. 



B. A general criterion for spinon confinement 

More importantly, however, the exact models we in- 
troduced above provide information about the models of 
SU(n) spin chains with representations corresponding to 
Young tableaux with a number of boxes divisible by n, 
i.e., models for which the Affleck-Licb theorem is not 
applicable. We have studied two SU(3) models belong- 
ing to this family in Sec. VIII, the rep. 10 VBS and the 
rep. 8 VBS, and found that both have confined spinons 
or colorons and hence display a Haldane-type gap in the 
spectrum. 

In this section, we will argue that models of antifer- 
romagnetic chains of SU(n) spins transforming under a 
representation corresponding to a Young tableau consist- 
ing of a number of boxes A divisible by n generally possess 



a Haldane-type gap due to spinon confinement forces. 

We should caution immediately that our argument is 
based on several assumptions, which we consider reason- 
able, but which we are ultimately unable to prove. 

The first, and also the most crucial, is the assump- 
tion that the question of whether a given model supports 
free spinon excitations can be resolved through study 
of the corresponding VBS state. This assumption def- 
initely holds for SU(2) spin chains, where the MG model 
for S = ^ indicates that the spinons are free, while the 
AKLT model for S = 1 serves as a paradigm for spinon 
confinement and hence the Haldane gap. The general 
conclusions we derived from our studies of the SU(3) 
VBSs above rely on this assumption. The numerical re- 
sults we reported on the rep. 8 VBS provide evidence 
that this assumption holds at least for this model. 

Let us consider an SU(n) spin chain with spins trans- 
forming under a representation corresponding to Young 
tableau consisting of L columns with Ki < K2 < . ■ . < 
Kl < n boxes each, 



^ I 



(81) 



with a total number of boxes 

L 

divisible by n. We denote the Dynkin coordinates of this 
representation by (/ii, ^2, • • ■ , /in-i), which implies 

n-l 



Note that this representation is, by definition, given by 
the maximally symmetric component of the tensor prod- 
uct of the individual columns. 



[k1, K2, . . . , Kl] = S{[ki] (g) [K2] 



])■ (82) 



For convenience, we denote the (") dimensional repre- 
sentation [k] in this subsection as ki = [k/] = (T). 

Since A is divisible by n, it will always be possible to 
obtain a singlet from the complete sequence of represen- 
tations Kl, K2, . . . , Kl by combining them antisymmet- 
rically. To be precise, when we write that we combine 
representations Ki and K2 antisymmetrically, we mean 
we obtain a new representation [ki -|- K2] by stacking the 
two columns with ki and K2 boxes on top of each other 
if Kl + K2 < n, and a new representation [ki + K2 — n] if 
Kl + K2 > n. In equations, we write this as 



A{[ki] «) [K2]) 



[ki + K2] for Ki + K2 <n ^gg^ 

[ki + K2 — 7l] for Kl + K2 > n 
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Following the notation used above, we indicate the an- 
tisymmetric combination of representations ni by a line 
connecting them. In particular, we depict the singlet 
formed by combining Ki, K2, ■ ■ ■ , k,l on different lattice 
sites as 



(84) 



The understanding here is that we combine them in the 
order indicated by the line, i.e., in (84) we first combine 
Ki and /«2, then we combine the result with K3, and so 
on. Depending on the order of the representations ft; on 
the line, we obtain different, but not necessarily orthog- 
onal, singlets. We assume that it is irrelevant whether 
we combine the representations starting from the left or 
from the right of the line, as the resulting state will not 
depend on it. 

In general, it will be possible to construct a number 
< \/n oi singlets out of various combinations of the k's, 
one for each block of k's for which the values of k add up 
to n as we combine the representation in the order de- 
scribed above. In this case, we will be able to construct 
one VBS for each singlet, and subsequently combine them 
at each site symmetrically to obtain the desired repre- 
sentation (81). The argument for spinon confinement we 
construct below will hold for each of the individual VBSs, 
and hence for the combined VBS as well. It is hence suf- 
ficient for our purposes to consider situations where the 
entire sequence Ki, K2, . . . , kl is needed to construct a 
singlet. 

A possible VBS "ground state" for a representation 
corresponding to a Young tableau with L = 4 columns is 
depicted below. 

one site 




(85) 



projection onto representation 

K2, K3, K4] 

In general, there are as many inequivalent VBS "ground 
states" as there are inequivalent ways to order the repre- 
sentations K.2, ■ ■ ■ , i.e., the number of inequivalent 
VBS "ground states" is given by 

L! 



Ail! • A*2! • 



,1' 



To give an example, the following VBS 
one site 




is inequivalent to the one above if and only if ki 7^ K2. 
Note that all these VBS "ground states" states are trans- 
lationally invariant. We expect some of these states, but 
not all of them, degenerate in energy for the appropri- 
ate parent Hamiltonian, and have accordingly written 
"ground states" in quotation marks. For example, if 
we form a SU(4) VBS with representation [ki,K2,'*3] = 
[1,1,2], the combination 



might yield a state with a lower energy for the appropri- 
ate Hamiltonian than the state formed by combining 



Simple examples where we have only two inequivalent 
VBS ground states are provided by the matrix product 
states discussed in Sees. VIII C and IXC. 

We will now argue that the elementary excitations of 
the corresponding VBS models are always confined. To 
this end, we first observe that any domain wall between 
translationally invariant "ground states" consists of a to- 
tal of m • n representations ki (m integer). To illustrate 
this, consider a domain wall between the "ground states" 
depicted in (86) and (85): 




projection onto representation 

K2, K3, K4] 



(86) 



In the example, the domain wall consists of reps. k,2, 
K2, /«3, and /«4. If the translationally invariant states 
on both sites are true ground states, the domain wall 
is likely to correspond to two elementary excitations: a 
rep. hi spinon consisting of an antisymmetric combina- 
tion of a /«2, a /^s, and indicated by the line 
in the drawing, and another rep. K2 spinon. The reason 
we assume that K2, /«3, and K4 form a rep. k,i is simply 
that this combination is present in both ground states 
on either side, and hence bound to be the energetically 
most favorable combination. The second K2, however, 
is not part of this elementary excitation, as combining 
it antisymmetrically with the others (i.e., the ki) would 
induce correlations which are not present in the ground 
state. We hence conclude that the second K2 is an el- 
ementary excitation as well. The domain wall depicted 
above consists of a spinon transforming under rep. hi 
and an anti-spinon transforming under rep. K2. 

The next step in our argument is to notice that the 
spinon and the complementary particle created simul- 
taneously which may either be its anti-particle or some 
other spinon, are confined through a linear potential. To 
see this, we pull them apart and look at the state inbe- 
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tween (color online): 



to the confining potential in (87) above (color online): 



K2-spinon 



Ki-spinon 



Ki-spmon 



Ki-spinon 




energy cost oc distance 

(87) 

The state between spinon and anti-spinon is not trans- 
lationally invariant. In the example, the unit cell of this 
state is depicted in red and consists of two regular bonds 
with three "singlet lines" between the sites, one stronger 
bond with four lines (which cross in the cartoon), and one 
weaker bond with only two lines. If we assume that the 
two states (86) and (85) on both sides are true ground 
states, it is clear that the irregularities in the strength 
of the bond correlations will cause the state between the 
spinon and anti-spinon to have a higher energy than ei- 
ther of them. This additional energy cost will induce 
a linear confinement potential between the spinous, and 
hence a linear oscillator potential for the relative motion 
of the particles. As in the models studied above, the 
Haldane gap corresponds to the zero-point energy of this 
oscillator. 

If one of the "ground states" to the left or to the right of 
the domain wall is not a true ground state, but a transla- 
tionally invariant state corresponding to a higher energy 
than the ground state, there will be a confining force be- 
tween this domain wall and another domain wall which 
transforms the intermediate "ground state" with a higher 
energy back into a true ground state. This force will be 
sufficient to account for a Haldane gap, regardless of the 
strength or existence of a confinement force between the 
constituent particles of each domain wall. 

The lowest-lying excitations of a representation 
K2, . . . , kl\ spin chain, however, will in general not 
be domain walls, but spinous created by breaking up one 
of the singlets (84) in a ground state. This is illustrated 
below for the ground state (85): 




In the example, we have created a spinon transforming 
under rep. Ki and its anti-particle, a spinon transform- 
ing under rep. Hi. This is, however, irrelevant to the 
argument — we may break the singlet in any way we like. 
The important feature is that we obtain, by construction, 
at least two excitations, and that these are confined. In 
our specific example, the confining potential is equivalent 




energy cost oc distance 

We leave it to the reader to convince him- or herself that 
the conclusions regarding confinement drawn from the 
simple examples studied here hold in general. 



C. Models with confinement through interactions 
extending beyond nearest neighbors 

Let us briefly summarize the results obtained. The 
SU(n) models we have studied so far fall into two cat- 
egories. The models belonging to the first — the trimer 
chain, the 6 VBS, and the n-mer chain — have n degener- 
ate ground states, which break translational invariance 
up to translations by n lattice spacings. The Young 
tableaux describing the representations of SU(n) at each 
site consist of a number of boxes A which is smaller than 
n, and hence obviously not divisible by n. The models 
support deconfined spinon excitations, and hence do not 
possess a Haldane gap in the spectrum. The Hamiltoni- 
ans of these models require interactions between n + \ 
neighboring sites along the chain. Even though the 
Afflcck-Lieb theorem is not directly applicable to the 
models we constructed above, it is applicable to SU(n) 
spin chains with spins transforming under the same rep- 
resentations. Like the VBS models, the theorem suggests 
that there is no Haldane gap in this family of models. 

The models belonging to the second category — the 10 
VBS, the 8 VBS, the representation (n,0, ...,0) VBS, 
and the m-rn matrix product state — have translationally 
invariant ground states. The ground states are unique 
for some models, like the 10 and the (n, 0, . . . , 0) VBS, 
and degenerate for others, like the 8 VBS. The Young 
tableaux describing the representations of SU(n) at each 
site consist of a number of boxes A which is divisible by 
n. The Affleck-Lieb theorem is not applicable to models 
of this category. The spinon excitations for this cate- 
gory of models are subject to confinement forces, which 
give rise to a Haldane gap. The parent Hamiltonians 
for these models require interactions between nearest- 
neighbor sites only. 

At first glance, this classification might appear com- 
plete. Further possibilities arise, however, in SU(7i) spin 
chains where number of boxes A the Young tableau con- 
sists of and n have a common divisor different from n, 
which obviously requires that n is not a prime number. 
In this case, it is possible to construct VBS models in 
which the ground state breaks translational invariance 
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only up to translations by n/ q lattice spacings, where q 
denotes the largest eommon divisor of A and n such that 
q < n. This implies that the ground state of the ap- 
propriate, translationally invariant Hamiltonian will be 
n/q-fold degenerate. In the examples we will elaborate 
on below, the parent Hamiltonians for these models re- 
quire interactions between -I- 1 sites, a feature we con- 
jecture to hold in general. The spinon excitations of these 
models are confined, even though the AfHeck-Lieb the- 
orem states that the nearest-neighbor Heisenberg chain 
of SU(n) spins transforming under these representations 
is gapless. (We implicitly assume here that the ground 
states of the SU(n) nearest-neighbor Heisenberg chains 
are non-degenerate.) Let us illustrate the general fea- 
tures of this third category of models with a few simple 
examples. 

(1) Consider an SU(4) chain with spins transforming 
under the 10-dimensional representation 

(2,0,0) = m . 

Following the construction principle of the 6 VBS of 
SU(3), we find that the two degenerate VBS states il- 
lustrated through 



one site 




T 

projection onto rep. (2,0,0) (88) 
are exact zero-energy ground states of 

i=l ^ ^ 

(89) 

which contains next-nearest neighbor interactions. 

The example illustrates the general rule stated above. 
The largest common divisor of n = 4 and A = 2 is 9 = 2. 
This implies ^ = 2 and hence two degenerate VBS states 
which break translational invariance up to translations by 
two lattice spacings. The parent Hamiltonian requires 
interaction between three neighboring sites. 

According to the AfHeck-Lieb theorem, models of an- 
tiferromagnetic SU(4) chains of representation (2,0,0) 
with nearest-neighbor Heisenberg interactions and non- 
degenerate ground states are gapless in the thermody- 
namic limit, which implies that the spinous are decon- 
fined. In all the models we have studied in previous sec- 
tions, the conclusions drawn from the Affleck-Lieb the- 
orem were consistent with those drawn from our exact 
models. For the present model, however, they are not 
consistent. 

Specifically, we strongly conjecture that the spinous in 
the (2,0,0) VBS are confined. This conjecture is based 
on the reasonable assumption that the elementary exci- 
tations of the model transform as either the fundamental 
representation 4 = (1, 0, 0) of SU(4) or its conjugate rep- 
resentation 4 — (0,0,1). This assumption implies that 



a single domain wall in one of the 4-mer chains used to 
construct the VBS state (88) shifts this chain by one lat- 
tice spacing. If we assume a ground state to the left of 
the spinon, the state on to the right will have a higher 
energy for the next-nearest neighbor Hamiltonian (89), 
as illustrated below. 



spinon spinon 
4 4 




energy cost oc distance (90) 

To recover the ground state, a second domain wall is 
required nearby, which is bound to the first by a linear 
potential. 

Our conclusion is not in contradiction with the rigor- 
ous result of Affleck and Lieb, as they explicitly restrict 
themselves to models with nearest-neighbor interactions. 
If we had only nearest-neighbor interactions, the energy 
expectation value in the region between the domain walls 
would not be higher than in the ground state, as one can 
see easily from the cartoon above — the sequence of al- 
ternating links would merely shift from (strong, weak, 
strong, weak) to (strong, strong, weak, weak). 

(2) The situation is similar for an SU(6) chain with 
spins transforming under the 56-dimensional representa- 
tion 

(3,0,0,0,0) = nJ3 ■ 

With n = 6 and A = 3, we have again ^ — 2. Accord- 
ingly, we find that the two VBS states illustrated through 



one site 




projection onto rep. (3,0,0,0,0) (91) 

are exact ground states of a parent Hamiltonian contain- 
ing up to next-nearest-neighbor interactions only, and 
that the spinon excitations are confined. 

(3) As a third example, consider an SU(6) spin chain 
with spins transforming under the 21-dimensional repre- 
sentation 

(2,0,0,0,0) ^ m ■ 

This imphes f = 3. We find that the three VBS states 
illustrated by 



one site 




projection onto rep. (2, 0, 0, 0, 0) (92) 

are exact ground states of a parent Hamiltonian involving 
the quadratic Casimir of total spin of four neighboring 
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sites, 
with 



N 



(2,0,0,0,0) VBS 



(93) 



r(4) 



(4)' 



44 



(4) 



50 

y 



(94) 

These VBS states break translational symmetry only up 
to translations by three lattice spacings. The spinons of 
this model are again confined through a linear potential, 
as illustrated below. 



spinon 
6 



spmon 
6 

o 
o 



energy cost oc distance 



(95) 



The conclusions wc have drawn for this VBS model rest 
on the assumption that the quadratic Casimirs of the 
representations contained in the tensor product shown in 
the dashed box in (92) as well as in the tensor product 
one obtains if one shifts this box by one lattice spacing 
to the left or to the right are smaller than the largest 
Casimir contained in the tensor product shown in the 
dotted box in (95). We have verified the vahdity of this 
assumption for the (2.0,0,0,0) VBS model we considered 
here, but not shown it to hold for similar models with 
larger n or X. 

(4) Finally, consider an SU(6) spin chain with spins 
transforming under the 70-dimcnsional representation 

(1,1,0,0,0) = p . 

Thus we have once again ^ = 2. In a notation similar to 
the one introduced for the 8 VBS, 

□ =o, g =o, 

the two degenerate VBSs are illustrated by 
one site 




projection onto rep. (1, 1, 0, 0, 0) 



(96) 



are exact ground states of a parent Hamiltonian involving 
the quadratic Casimir of the total spin of three neighbor- 
ing sites. 



N 



1,0,0,0) VBS — Hi 



(97) 



with 



= ( [ji'^Y- 20] ( (J^^- 70] ( fjf V- 540 



(98) 



The states (96) are not the only VBSs one can form. 
Other possibilities like 



or 




however, contain additional representations for the to- 
tal SU(6) spin of three neighboring sites, and are hence 
expected to possess higher energies. 

Spinon excitations transforming under the 6- 
dimensional rep. (0,0,0,0,1) arc linearly confined 
to spinons transforming under the 15-dimensional rep. 
(0,0,0,1,0): 



spinon 
(0,0,0,0,1 



spmon 
(0,0,0,1,0) 




energy cost oc distance 



(99) 



The VBS configuration we have drawn between the two 
spinons in (99) constitutes just one of several possibilities. 
We expect, however, that this possibility corresponds to 
the lowest energy among them for the Hamiltonian (97). 
This concludes our list of examples. 

The models introduced in this subsection are interest- 
ing in that they provide us with examples where spinon 
confinement, and with the confinement the existence of 
a Haldane gap, is caused by interactions extending be- 
yond nearest neighbors. The conclusion drawn from the 
Affleck-Lieb theorem for SU(n) models with spins trans- 
forming under representations we have labeled here as 
the "third category" hence appear to hold for models 
with nearest-neighbor interactions only, to which the the- 
orem is applicable. For these models, the theorem states 
that the spectrum is gapless, which according to our un- 
derstanding implies that the models support deconfined 
spinon excitations. The examples wc have studied, by 
contrast, suggest that models with longer-ranged inter- 
actions belonging to this category exhibit confinement 
forces between the spinon excitations and hence possess 
a Haldane gap. 

It is worth noting that even though in the examples 
we elaborated here A < n, we expect our conclusions to 
hold for models with A > n as well. To sec this, let 
m > be an integer such that nm < A < n{m +1). We 
can now construct a first VBS with spinon confinement 
using nm boxes of the Young tableau and combine it 
with a second by projection on each side with a second 
VBS constructed from the remaining A' = A — nm boxes. 
Since the spinons of the first VBS are always confined 
and hence gapped, the final VBS will support deconfined 
spinons if and only if the second VBS will support them, 
which in turn will depend on the largest common divisor 
q' of A' and n as well as the range of the interaction. 



21 



Since the largest common divisor g of n and A is equal to 
g', there is no need to think in terms of A' and g'. The 
conclusions regarding confinement and the Haldanc gap 
will not depend on the distinction between A and A'. 

D. The different categories of models 

In this section, we used the rules emerging from the 
numerous examples we studied to argue that models of 
SU(ri) spin chains in general fall into three categories. 
The classification depends on the number of boxes A the 
Young tableau corresponding to the representation of the 
individual spins consists off, as follows: 

1. If A and n have no common divisor, the models will 
support free spin excitations and hence not exhibit 
a Haldane gap. 

The general reasoning here is simply that the VBS 
states in this category break translational invari- 
ancc up to translations by n lattice spacings, and 
that there are (at least) n degenerate VBS ground 
states to each model. Spinous transforming un- 
der representations of Young tableaux with an ar- 
bitrary number of boxes can be accommodated in 
domain walls between these different ground states. 
Consequently, the spinous are deconfined. 

2. If A is divisible by n, the general argument we have 
constructed in Sec. XB indicates that the model 
will exhibit spinon confinement and hence a Hal- 
dane gap. 

3. If A and n have a common divisor different from 
71, the examples studied in Sec. XC suggest that 
the question of whether the spinous are confined or 
not depends on the details of the interactions. If 
q is the largest common divisor of A and n, inter- 
actions ranging to the -^-th neighbor were required 
for spinon confinement in the models we studied. 
The Affleck-Lieb theorem [63], on the other hand, 
tells us that SU(n) chains with nearest-neighbor 
Heisenberg interactions belonging to this category 
are gapless if the ground states are non-degenerate. 

Note that the second category is really just the special 
case g = n of the third: with ^ = 1, nearest- neighbor 
interactions are already sufficient for spinon confinement 
and a Haldane gap. 

The conclusion that interactions ranging to the -^-th 
neighbor are required for spinon confinement and a Hal- 
dane gap, however, is not universally valid. A counter- 
example is provided by the extended VBSs (XVBSs) in- 
troduced by Affleck, Arovas, Marston, and Rabson [75]. 
In this model, each site effectively takes the role of 
two neighboring sites when a VBS is constructed, and 
nearest-neighbor interactions are already sufficient to 
cause spinon confinement. We briefly review this model 
in the following subsection. 



E. A counter-example: the SU(4) representation 6 
extended VBS 

In an article devoted to quantum antiferromagnets 
with spins transforming under the self-conjugate repre- 
sentations of SU(2n), Affleck et al. [75] introduced an 
extended VBS for the six-dimensional SU(4) represen- 
tation (0,1,0). This representation is, with regard to 
the number of boxes the corresponding Young tableau 
consists of, not distinguishable from the symmetric rep- 
resentation 10 considered in Sec. XC, as both tableaux 
consist of two boxes: 

□ ® □ = m © 

4 4 10 6 (100) 

Since the two boxes are combined antisymmetrically for 
rep. 6, a VBS constructed along the lines of Sec. XC 
would no longer provide a paradigm for antiferromag- 
netic spin chains of the corresponding representation in 
general. Affleck et al. [75] have constructed an extended 
VBS, which is illustrated in the following cartoon: 

Here each small circle represents a fundamental represen- 
tation 4 of SU(4) (a box in the Young tableau), and each 
large circle a lattice site. The lines connecting four dots 
indicate that these four fundamental representations are 
combined into an SU(4) singlet. The total spin of two 
neighboring sites in this state may assume the represen- 
tations 

4® 4 =10 15, (102) 

while combining two reps. 6 on neighboring sites in gen- 
eral yields 

6(8)6 = 10 15 020. (103) 

To construct a parent Hamiltonian for the XVBS (101), 
it is hence sufficient to sum over projectors onto rep. 20 
on all pairs of neighboring sites. Using our conventions 
(Affleck et al. [75] have normalized the eigenvalues of the 
quadratic Casimir operator to C^y^^^ (adjoint rep.) = 8), 
the parent Hamiltonian takes the form 

" / 1 5 \ 

i^xvBS = + 3 i-J'^J^+if + ) . (104) 

where the operators Jf, a = 1, . . . , 15, are given by 6 x 6 
matrices. Note that the ground state is two fold degen- 
erate, as it breaks translational symmetry modulo trans- 
lations by two lattice spacings. 

We conjecture that the lowest lying excitation is a 
bound state consisting of two fundamental reps 4, which 
most likely are combined antisymmetrically into a rep. 6: 

^^^^^^^^^^^^ 
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The SU(4) rep. 6 XVBS provides us with an example 
where a nearest-neighbor Hamiltonian is sufBcient to in- 
duce spinon confinement and a Haldane gap, even though 
the largest common divisor of n = 4 and A = 2 is g = 2, 
i.e., interactions including = 2 neighbors would be re- 
quired following the rules derived from the examples in 
Sec. XC. The reason for this discrepancy is that in the 
XVBS model considered here, each site effectively takes 
the role of two neighboring sites. Affleck et al. [75, 76] 
conjecture that the ground state of the SU(4) rep. 6 
nearest-neighbor Hcisenbcrg model is, like the XVBS re- 
viewed here, two-fold degenerate, which implies that the 
Affleck-Lieb theorem is not applicable. 

This example is valuable in showing that it is advisable 
to explicitly construct the VBS for a given representation 
of SU(n) in order to verify the applicability of the general 
rules motivated above. 



XI. CONCLUSION 

In the first part of this article, we have formulated 
several exact models of SU(3) spin chains. We intro- 
duced a trimer model and presented evidence that the 
elementary excitations of the model transform under the 
SU(3) representations conjugate to the representation of 
the original spin on the chain. We then introduced three 
SU(3) valence bond solid chains with spins transforming 
under representations 6, 10, and 8, respectively. We ar- 
gued that of these four models, the coloron excitations 
are confined only in the 10 and the 8 VBS models, and 
that only those models with confined spinous exhibit a 
Haldane gap. We subsequently generalized three of our 
models to SU(7i), and investigated again which models 
exhibit spinon confinement. 

Finally, we used the rules emerging from the numer- 
ous examples we studied to argue that models of SU(n) 
spin chains in general fall into three categories with re- 
gard to spinon confinement and the Haldane gap. These 
are summarized in the previous subsection. The results 
rely crucially on the assumption that the conclusions we 
obtained for the VBS models we studied are of general 
validity. This assumption certainly holds for the corre- 
sponding SU(2) models, and appears reasonable on phys- 
ical grounds. Ultimately, however, it is only an assump- 
tion, or at best a hypotheses. 

On a broader perspective, we believe that the models 
we have studied provide further indication that SU(n) 
spin chains are an equally rich and rewarding subject of 
study as SU(2) spin chains have been since Bethe. 
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APPENDIX A: A PROPOSAL FOR AN 
EXPERIMENTAL REALIZATION OF SU(3) SPIN 
CHAINS IN AN OPTICAL LATTICE 

In this appendix, we wish to describe a proposal for 
an experimental realization of SU(3) spin chains. The 
most eligible candidate for such experiments are ultra- 
cold gases in optical lattices. Recently, these systems 
have become an interesting playground for the realiza- 
tion of various problems of condensed matter physics, 
such as the phase transition from a superfiuid to a Mott 
insulator [77, 78], the fcrmionic Hubbard model [79-81], 
and SU(2) spin chains [82, 83]. In particular, the Hamil- 
tonians for spin lattice models may be engineered with 
polar molecules stored in optical lattices, where the spin 
is represented by a single-valence electron of a heteronu- 
clear molecule [84, 85]. 

In a most naive approach, one might expect to realize 
an SU(3) spin (at a site in an optical lattice) by using 
atoms with three internal states, like an atom with spin 
S" = 1. If we now were to interpret the S*^ = +1 state 
as SU(3) spin "blue", the 5"^ = state as "red", and the 
S"^ = —1 state as "green", however, the SU(3) spin would 
not be conserved. The SU(2) algebra would allow for 
the process [ + 1,— 1) |0, 0), which in SU(3) language 
corresponds to the forbidden process [b, g) |r, r)- 

A more sophisticated approach is hence required. One 
way to obtain a system with three internal states in 
which the number of particles in each state (i.e., of 
each color) is conserved is to manipulate an atomic sys- 
tem with total angular momentum = 3/2 (where 
F = Sci + iorb + S'nuc iucludes the internal spin of the 
electrons, the orbital angular momentum, and the spin of 
the nucleus) to simulate an SU(3) spin. The important 
feature here is that the atoms have four internal states, 
corresponding toi^^ = — |,— ^,-|-^,-|-|. For such atoms, 
one has to suppress the occupation of one of the "middle" 
states, say the — state, by effectively lifting it to 
a higher energy while keeping the other states approxi- 
mately degenerate. This can be accomplished through a 
combination of an external magnetic field and two care- 
fully tuned lasers, which effectively push down the en- 
ergies of the i^^ = — I and the = +i states by 
coupling these states to states of (say) the energetically 
higher F = 5/2 multiplet (see Fig. 15). At sufficiently 
low temperatures, we are hence left with a system with 
three internal states F^ = — which we may 
identify with the colors "blue" , "red" , and "green" of an 
SU(3) spin. In leading order, the number of particles of 
each color is now conserved, as required by SU(3) sym- 
metry. For example, conservation of F^ forbids processes 
in which a "blue" and a "green" particle turn into two 
"red" ones, |b, g) |r,r). Higher order processes of the 
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"blue" "red" "green" 

FIG. 15: Effective lifting of the i^' = - i state. 



kind |b,g, g) kjrjr) are still possible, but negligible 
if the experiment is conducted at sufficiently short time 
scales. 

If one places fermionic atoms with an artificial SU(3) 
spin engineered along the lines of this or a related pro- 
posal in an optical lattice and allows for a weak hopping 
of the atoms on the lattice, one has developed an exper- 
imental realization of an SU(3) Hubbard model. If the 
energy cost U of having two atoms on the same lattice 
site is significantly larger than the hopping t, and the 
density is one atom per site, the system will effectively 
constitute an SU(3) antiferromagnet. The dimension of 
this antiferromagnet will depend on the optical lattice, 
which can be one-, two-, or three-dimensional. 

In principle, the above proposal can be generalized to 
SU(n), even though the experimental obstacles are likely 
to grow "exponentially" with n. Besides, it is far from 
clear that such an endeavor is worthwhile, as all the non- 
trivial properties of SU(n) are already present in SU(3) 
chains (while SU(2) constitutes a special case). 



Explicitly, the non-vanishing structure constants are 

given by = i, /"7 = /246 ^ J257 ^ J345 ^ _jl56 ^ 

-/367 = j/2, = /678 = ^^/3/2, and 45 others 

obtained by permutations of the indices. 



APPENDIX C: EIGENVALUES OF THE 
QUADRATIC CASIMIR OPERATOR 

The eigenvalues of the quadratic Casimir operator for 
representations Cgu^„j(/ii, ^2, ■ • ■ , Mn-i) of SU(7t.) up to 
n = 6 are given by: 

Ciu(2)(/^) = i(M^ + 2M) = f + 1) 

'^su(3)(/^i'/^2) = K^ii + MiA*2 + M2 + 3^1 + 3^2) 

^su(4)(/^i'M2,Ai3) |(3Ai?+4^i + 3/.ti + 4^iAi2 
+2fiifi3 + 4^2/i3 + 12mi + 16m2 + 12^3) 

<^su(5) (mi ' M2 , M3 , M4) = I (2Aii + 3pi| -I- S^ij + 2^1 

+ 3Aii/^2 + 4^2/^3 + i^J■3^J"l + 2MiAi3 + ^J'l^J■4 

+2^2fJ.4 + -I- I5fj.2 + 15/^3 + IO//4) 

'"SU(6)(/^1' ^2,^3:^4,^*5) = 

^ {5^il + 81^4 + 9/^3 + 8^*1 + 5^^ 

+ 8^i^2 + 12Ai2M3 + 12/i3Ai4 + 8^4M5 
+ 4/ii^4 + 6^i/i3 + 8^2/"4 + 6/i3^5 

+ 4/i2^5 + 2^i/i5 -I- 30/.J1 + 48^2 
+Mi^L3 + 48^*4 + 30^5) 



APPENDIX B: GELL-MANN MATRICES 

The Gell-Mann matrices are given by [67, 
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They are normalized as tr (A''A'') = 26ab and satisfy the 
commutation relations [A°, A''] = 2/"''^A'^. The structure 
constants /"''^ are totally antisymmetric and obey Ja- 
cobi's identity 



j^abc jcde jbdc jcae 



0. 



The general method to obtain these and further eigen- 
values for n > 6 requires a discussion of representation 
theory [74] at a level which is beyond the scope of this 
article. 

The dimensionality of a representation 
(/ii, /j,2, . . . , /i„_i) is determined by the so-called 
Hook formula [67] 



dim 



(Cl) 



where Xi = Mi for i = 1, . . . ,n. In particular, it 

yields for n = 2, 3, 4: 

dinisu(2)(M) = + 1 

dimsu(3)(Mi,/"2) = ^{f^i + l)(/i2 + l)(/^i + M2 + 2) 

dhTlsu(4)(Ml,M2,M3) = + 1)(M2 + 1)(M3 + 1) 

i^J■l + ^J■2+ 2)(/i2 + Af3 + 2)(/il + M2 + Ai3 + 3) 
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